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Abstract 

We use chiral perturbation theory to evaluate the scattering amphtude for the process 
'K~^K~ TT'^K~ at leading and next-to-leading orders in the chiral counting and in the 
presence of isospin breaking effects. We also discuss the influence of the latter on the 
combination of the 5- wave ttK scattering lengths which is relevant for the 2S — 2P energy 
levels shift of Ktt atoms. 
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1 Introduction 



The study of hadronic atoms has become a very active field. Many experiments are devoted to 
measure the characteristics of such atoms with high precision [|I|, ^, ^, |^ . These experimental 
results carry a theoretical interest since they permit a direct access to hadronic scattering 
lengths providing, in this way, valuable informations concerning the fundamental properties of 
QCD at low energy. For instance, the presently running DIRAC experiment aims at measuring 
the pionium lifetime r with 10% accuracy This would allow one to determine the difference 



with 5% precision by means of the Deser-type [[ 

ao) 



relation 



oc a, 



(1.1) 

where aj denotes the /-wave nn scattering length in the channel with total isospin I. On the 
other hand, chiral perturbation theory (ChPT) H, P, predictions for the scattering lengths 
have reached a precision amounting to 2% Once the final results from DIRAC are available, 
ChPT will therefore be subjected to a serious test. Before confronting the experimental deter- 
mination to the ChPT prediction, it is desirable to get all sources of corrections to the relation 
(|1.1| ), valid in the absence of isospin breaking, under control. In this direction, bound states 
calculations were performed using different approaches, like potential scattering theory |1^, p!3[1 , 
3D-constraint field theory [14], Bethe-Salpeter equation |15| and non-relativistic effective la- 



grangians |T^. For a review on the subject and a comparison between the various methods 



we refer the reader to Within the framework of non-relativistic effective lagrangians, the 
correct expression of relation ( |1 . 1|) which include all isospin breaking effects at leading order 
(LO) and next-to-leading order (NLO) was found in |jl9[ to be 



r 



-1 



4M!o 



aY'A^l + K). 



;i.2) 



In the preceding equation, a stands for the fine-structure constant, A and K possess the 



following expansions 
A = 



in powers of the isospin breaking parameter k G [a, (m^ 



;00 



1.3) 
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K = U^-l){a', + 2alY-^i\na-l){2aUal)+oin). (1.4) 
The quantity of interest, 

- tJ- Re v4^^7'°° = '^0 - «o + hi{md - rriuf + , (1.5) 

represents the real part of the ti^-k^ -k^tx^ scattering amphtude at order k, calculated at 
threshold within ChPT to any chiral order and from which we subtract the singular pieces 
behaving like and Ing whith q being the centre-of-mass three-momentum. While hi van- 
ishes, the coefficient h2 was calculated in at (9(e^p^) where p stands for a typical external 
momentum and e for the electric charge. 

In the first DIRAC proposal |[l|], it was also planned to measure the pionium 25 — 2P energy 
levels shift Ai?. The possibility to perform such a measurement was discussed in [^. A 



simultaneous measurement of r and AE would allow to pin down Oq and separately, since p2| 

oc 2a|] + ag. (1.6) 



Bound states calculation of the isospin breaking corrections to ( |1.6| ) were done in ||T3| using 
potential scattering theory. Non-relativistic effective lagrangian calculations concerning Ai? 
are not available. One might however expect that they will involve the quantity Re 
for the corresponding vr+vr^ tt+tt^ process. Electromagnetic corrections to the scattering 
amplitude for the process in question have been obtained at O(e^p^) in p3| . 



The proposal [|l[ was followed by another one extending the DIRAC project by considering 
the possibility of measuring the characteristics of Ktx atoms with 20% accuracy |^. The 
determination of the lifetime and energy levels shift for Ktx atoms will give access to the S'-wave 

1/2 3/2 

TxK scattering lengths Oq and Oq allowing to test ChPT in the three-flavor sector. Likewise, 
it is to be expected that the relevant quantities in the flnal expressions for the lifetime and 
the levels shift will involve the scattering amplitudes for tx'K'^ —>■ n^K^ and n'^K'' —>■ tx^K~ 
respectively. In the case of the former, isospin breaking corrections have already been discussed 
in [p4|, p5|]. The aim of the present work is to provide a similar treatement for the latter. 



The paper is organized as follows. In Sec. |^, the scattering amplitude for the process ti^K^ — > 
Ti^K^ will be calculated at NLO including isospin breaking effects and ignoring the emission 
of real soft photons. Being lengthy, the expression for the scattering amplitude is displayed in 
App. ^. We pursue with Sec. ^ where the analytic expressions for and Og^^ at NLO are 
given. We then discuss the sensitivity of the scattering lengths to the size of the next-to-next-to 
leading order (NNLO) by evaluating them using various inputs for the low-energy constants 
(LEC's). The threshold expansion of the process in question is performed in Sec. ^ where the 
effects of isospin breaking on the scattering lengths are evaluated. Finally, App. |B| collects the 
expressions for the loop functions needed in the calculation. 



2 Charged pion and kaon elastic scattering 

The elastic scattering process 

7r+(p+) + K-{k_) ^ 7r+(p'+) + K-{k'_) , (2.1) 
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is studied in terms of the Lorentz invariant Mandelstam variables 



verifying the on-shell relation s + t + u = 2(M^± + M|,±). These variables are related to the 
centre-of-mass three momentum q and scattering angle 6 by 



s 



t = -2g2(l - cos^) , 

u = (^^M2±+g2_^M^±+g2^ - 2g2(l + cos ^) . (2.2) 



Let Ai^ and A^++'++ denote the respective scattering amplitudes for the process (|2.1| ) and 
its crossed channel reaction tt~^K~^ —>■ -k^K^ . Then, in the isospin limit, defined by the vanishing 
of both the electric charge e and the up and down quarks masses difference {tjIu = rn^ , e = 0), 
the following relations hold 

M^+'++{s,t,u) = T^^\s,t,u), (2.3) 

with (J = 1/2,3/2) being the isospin amplitudes. The processes under consideration are 
related by s ^ u crossing which constrains the by 

2T'/\s, t, u) = 3T3/2(u, t, s) - T''l\s, t, u) . (2.4) 

This is no more valid when isospin breaking effects, generated by 5 = rrid — and a = / (47r) 
are switched on. In this case s <-> u crossing is expressed as 

M^-''^-{s,t,u) = A<++'++(n,t,s). (2.5) 

The effect of a non-zero value for 6 can be analyzed fully by means of the strong sector chiral 



lagrangian constructed in |]26[ . Treating isospin violation of electromagnetic origin requires the 
extension of ChPT in order to include virtual photons. This can be done by building operators 
in which photons figure as explicit dynamical degrees of freedom. The electromagnetic sector 
chiral lagrangian, founded upon the chiral counting scheme 0{e) = 0{p), was presented at 
NLO in and [^. 



We shall calculate the scattering amplitude ( ^.5[ ) at NLO including isospin breaking effects of 
both strong and electromagnetic origin. Consistency demands that all of the following chiral 
orders should be present; p^, e^, 5, p^, e^p^, e^, 5e^, 5^. From naive dimensional estimation, 
we believe that the last three orders are beyond the accuracy we ask for and hence will be 



ignored in the following. Furthermore instead of S, our results will be expressed in terms of p6 



4 nis — 171 

which measures the rate of isospin violation with respect to the violation of SU{3). Using 
Feynman graph techniques, the amplitude ( |2.5D can be represented at NLO by the one-particle- 
irreducible diagrams depicted in Fig. |^. 
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(d) (e) (f) 

Figure 1: The various types of Feynman diagrams encountered in the charged nK scattering 
to one-loop order and ignoring Oi^e^). The born-type diagram is represented by (a). Besides 
the contribution from the LEC's (the full square), it contains the tree contribution including 
the effects of bare masses and wave functions renormalization. Diagram (b) represents the 
tadpole-type part of the amplitude. The s-, t- and w-channel parts are schematized by diagram 
(c) and its crossed. The one-photon contribution to the amplitude follows from diagrams (d) 
and (e). We refer to diagram (f) as the form- factor- type where the full circle is made explicit 
in Fig. 1^. 




Figure 2: The electromagnetic vertex function of a charged meson to one-loop order. The 
full square takes into account the contribution from the LEC's just as the tree contribution 
including the effect of wave function renormalization. Are shown only diagrams of 0{ep^). 
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Although the evaluation of these diagrams is a simple exercise in quantum field theory, it 
involves different masses, yielding, in that manner, lengthy expressions. These are displayed 
in App. 1^ where we keep the contributions of the various diagrams separated. It is useful to 
mention, at this stage, that these expressions are scale independent but infrared divergent. Since 
only observables are infrared safe, these infrared divergencies should cancel in the expression 
for the cross section where not only virtual photons but also real soft photons have to be 
taken into account. The 0{a) soft photon cross section corresponding to an arbitrary matrix 
element was calculated in |^9[ . We cheked that by applying the general formula derived in that 
reference to the present case, the terms in In so obtained cancel with the ones coming from 
our expressions. 



3 Scattering lengths 



We are interested by the S'-wave nK scattering lengths pO|, |32|] . To this end, it is convenient 
to introduce the partial wave amplitudes t[ defined in the s-channel by 



s, cos ( 



327r^(2/ + l)P;(cos^)tf(s) 



(3.1) 



where / is the angular momentum of the nK system and the P;'s are the Legendre polynomials. 
Near threshold the partial wave amplitudes can be parametrized in terms of the scattering 
lengths a[ and slope parameters bj . In the normalization ( |3.1|) , the real part of the partial wave 
amplitudes reads 

Retiis) = g2^[af + 6fg2 + 0(g^)] . (3.2) 
At NLO in ChPT no analytic expressions for the scattering lengths exist in the literature. 



This is not the case in other approaches such as heavy-kaon ChPT where the expansion 
parameter is M^/Mk in addition to M^/(47rFo) or M7^/(47rFo) for ChPT^. Obviously, at any 
order in both expansions, a matching between the two approaches is feasible. Getting analytic 
expressions for and ajj^^ is a straightforward matter. Using the expression of the isospin 
amplitude T^^^ given in g, Eqs. (|23)-(|23) together with (^ and (^) lead the following 
NLO expressions for the scattering lengths 



327ra, 



1/2 



1 + ^^(M2±+M|.±)L 



FttFk 

1 



1 



1 



9M2±(11M|± - 5M2±) In 



+ 2M|:±(9Mi± - 55M2±) In — |i + (36M^± + 11M|;±M2± - 9M^±) In 



+ mM^[A{L\ + LI) + 2(L3 - 2LI) - L; + 2{2Ll + L;)] + 



Ml. - 



X 



^AU along this paper, M-j^ and Mk represent respectively the pion and kaon masses in the isospin limit, 
the eta mass, Fq the coupling of Goldston bosons to axial currents in the chiral limit. 
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M2 



327ra' 



3/2 



99M^± In - 2(67M^± - 8M^±) In 
M^±Mk± ( 



M 



+ (24M|± 



5M^ 



M2 
/i2 



1 



1^ 



4 



F-kFk 
1 



+ 



+ 



(9M- 



M2, 



9Mi 



(IIM 



5Mi 



)ln 



2M2± 



In- 



+ (36M^± + llMl^Ml^ - 9M^± 



M2 



1 172 + 3847r2i3(M^±, -M^±) 



+ 46 087r2[4(L^ + L^) + 2(L3 - 21^) - L; + 2(2L^ + L^)] + 



99M2± In 



M 



2(67M^± - 8M^±) In + (24M|.± 



M2 
In^ 

/i2 



(3.3) 



The preceding expressions were obtained by using, in the NLO contributions, the Gell— Mann- 
Okubo relation ([4.5|) in the isospin hmit [^. and Fx stand for the pion and kaon decay 
constants respectively. Their NLO expressions in the isospin limit can be found in ^6[. For 
their numerical values we will use = 92.4 MeV |3^ and Fx = 1.22F^ [37, 38|. The Lj's 
are the LEC's weighting O(p^) operators in the effective lagrangian of and whose values 



are collected in Tab. with various experimental determinations. Finally, for compactness, we 
introduced the function B whose expression is displayed in App. |B[ For historical reasons [Q, 
the scattering lengths were defined in terms of the charged pion and kaon masses {Mt^± = 
139.570 MeV, M/^± = 493.677 MeV). Furthermore, F^ was renormalized as F^^Fx instead of 
F^. Nevertheless, if one whishes to adopt the second choice for the renormalization of Fq, the 
isospin limit of Eq. ([4.7|) can be used. 





set I 


set II 


set III 




0.46 ±0.23 


0.53 ±0.25 


0.43 ±0.12 




1.49 ±0.23 


0.71 ±0.27 


0.73 ±0.12 




-3.18 ±0.85 


-2.72 ± 1.12 


-2.35 ±0.37 




0±0.5 


0±0.5 


0±0.5 




1.46 ±0.2 


0.91 ±0.15 


0.97±0.11 


WLl{M,) 


0±0.3 


0±0.3 


0±0.3 


103L^(M,) 


1.00 ±0.20 


0.62 ±0.20 


0.60 ±0.18 



Table 1: Values of the L['s obtained in ||3^ by using large N^. arguments and fitting data 



from ^l|] to ChPT predictions at NLO (set I) as also at NNLO (set II). Set III is the equivalent 
of set II with data coming from the preliminary analysis of the E865 experiment. 



We applied several checks to expressions ( |3.3| ). They are scale independent, the scale depen- 
dence /i of the chiral logarithms is compensated by the one governing the renormalization group 
equations PB[ of the running couplings Ll{^) = L^. Even so, all of our expressions will be eval- 

1) to the fourth order in 



uated at the scale /i = Mp = 770 MeV. Expanding expressions 
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powers of M.,^±/Mk± we recover the combinations of scattering lengths obtained within the 
framework of heavy-kaon ChPT ^3|. Using the same inputs as in |Q we found again their 



numerical estimates. These expressions are consistent analytically and numerically with the 
combinations ^aH"^ + and — a^J"^ evaluated in ||2^ 



Fq renormalization 


«o 


set I 


set II 


set III 




1/2 
^0 

3/2 
«0 


0.214 ±0.016 
-0.0557 ±0.0166 


0.202 ±0.018 
-0.0660 ±0.0185 


0.203 ±0.013 
-0.0644 ±0.0133 


-^0^ = F-kPk 


1/2 
«0 

3/2 
«0 


0.192 ±0.011 
-0.0613 ±0.0113 


0.177 ±0.012 
-0.0651 ±0.0125 


0.179 ±0.009 
-0.0644 ±0.0090 



Table 2: ChPT predictions for the S'-wave t^K scattering lengths at NLO with two possible 
renormalization choices for Fq. Set I, set II and set III are defined in Table |l|. 



We shall now update the numerical evaluation of |]3^ for the scattering lengths. All inputs 



have been given exception of the eta mass to which we assign the value = 547.30 MeV. We 
will evaluate Cq^^ and Qq^'^ using each of the sets defined in Tab. Q with both renormalization 
choices for F^. The interest behind this is the following. Since the difference between F^^ and 
Fk in the NLO pieces is of 0{p^), then any difference in the values for the scattering lengths 
due to the renormalization choice for Fq could be viewed as an indication for the size of the 
NNLO corrections. With respect to the different sets of Tab. |I], set I and set II were obtained 
by fitting respectively the same experimental data to one- and two-loop ChPT predictions. It 
follows that the possible variation in the values for the scattering lengths due to the use of either 
the two sets measures the influence of the NNLO on the NLO which, in a way, reflects the rate 
of convergence of the chiral expansion in powers of m^. While Table |^ collects our numerical 
results for the scattering lengths, we will concentrate on two specific combinations, 2al^'^ + a^^'^ 
and Oq^^ — Qq^'^. For the latter, the variations are about ~ 9% due to the renormalization choice 
for Fo and ~ 1% (~ 5%) due to the choice of set II instead of set I with F^ = F^ {Fg = F^Fk). 
With regard to the former, the renormalization choice for Fq induces about ~ 17% of variation 
versus ~ 10% due to the choice of the set for both renormalization choices. Notice that up to 
10% of accuracy, the combination altered neither by the renormalization choice 

for Fq nor by the choice between set I and set II. This result has already been noticed in 



1/2 3/2 

where it was concluded that the theoretical prediction for Oq — Og is so clean that if the 
combination in question is measured accurately, the validity of standard three-flavor ChPT can 
be confirmed or not. 



4 Isospin violation at NLO 

The scattering lengths are well-defined quantities in the absence of radiative corrections. When 
virtual photons effects have to be taken into account, things must be handled with care. To 
prevent any confusion, we will perform step by step the threshold expansion of the amplitude 
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p.5| ) with the help of Apps. |^ and ^. When expanded in powers of q, the real part of ( |2.5| ) 



can be put in an analogous form as for the vrvr scattering case pO 



Re M-'-(s,t,u) = ^^^^^^ + e=(^)+ReX-;*- + 0(,). (4.1) 

ReA<-^-(M.«) = ^^^^f + .^(^)+ReA<r,r* + 0(,). (4.2) 

The first two terms in the right-hand side of equation ( [4. ID should be absorbed in the static 
characteristics of Kn atoms [|^, ^ ^ by means of a bound states treatement. The term in 



is due to the Coulomb photon exchanged between the scattered particles in diagram (d) of 
Fig. m. As to the term where the dependence on the Mandelstam variables is kept explicit, it 
corresponds to diagram (f) of Fig. |l| at tree level and contains, besides the dependence on 9, a 
singular piece behaving like The q- and 6'-independent terms. Re A^^~'^~ Re A4^^'^~^, 
constitute the central object of the current work. Their isospin limit is nothing else than the 
threshold value of equation ( p.3|) . Then, in order to keep a coherent notation, they will be 
denoted in the following by 

Re -Mttr^^ = 327rao(+T; +T) , 

where 

«o(+-;+-) = ^ (2ay' + a;;/') + Aao(+-; +-) , 

ao(++;++) = + Aao(++; ++) , (4.3) 



and the NLO expressions (|3.3|) were used for aj''^ and Oq^^. Let us turn towards the calculation 
of the isospin breaking quantities Aoq. To begin with, we express the isospin limit masses M^r 
and Mk figuring at NLO in terms of the charged and neutral ones by means of 

^ M^o , 2M|: ^ Mi± + M|o - l{M'^± - Ml,) , (4.4) 



where 7 takes into account any deviation from Dashen's theorem for which 7 = 1. Although 
our results correspond to 7 = 1, we will use the value 7 = 1.84 |^ as an indicator for their 
sensitivity to the violation of Dashen theorem |4^. Furthermore, for simplification, the modified 



Cell— Mann-Okubo relation 

3M^2 2(M|^± + M|o) - (2M2± - M^o) , (4.5) 

will be used everywhere except in the arguments of the chiral logarithms where the eta mass is 
assigned to its physical value. Nevertheless, the effect, on the final results, of ignoring Eq. ( |4.5| ) 
will be also considered. Next, since by convention, the isospin limit is defined in terms of the 
charged pion and kaon masses, the neutral ones should be replaced according to 

Ml, ^ Ml^ - A, , Ml, ^ - A^ . 

The advantage of such a procedure is that, when expanding ao in powers of Att and Aj^, the 
zero-th order in the expansion is automatically defined in terms of the charged masses and 
reproduces the expression for the corresponding combination of scattering lengths as given by 
(|4.3|) and (|3.3|). Although sufficient for our purposes, the expansion to first order in A^ and A^^ 
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is somewhat tedious especially at the level of the loop functions. The corresponding expansions 
of the latter are collected in App. |B[ Having all this, the last step is achieved by replacing 
and Ak in the NLO terms by their LO expressions 



which permit the following decomposition 
327rAao(+-;H— ) = 



F^Fk V3 ' Zoe2 



327rAao(++;++) = -^^^+-^5. 



F^Fk 73 ' Zoe2 



(4.6) 



Once more, the renormalization choice for Fq was fixed to Fq = F-j^Fk with the possibility of 
renormalizing as Fq = F^ offered by the relation 



F^Fi 



K 



1287r2F2 
A 



+ 



-\l + -,{Ml^-Ml^)Ll 

IT K TT 

In ^ - 2Ml^ In ^ - (4Mi± - Ml^) In 

1 / Ml M^^^ 

2L^--r— . (4 + 51n^-ln^-ln^ 

/i"' /i^ /i"' 

l + lnM + 21n^S^ 



Ak 
F3 



2LI 



1287r2 
1 



1287r^ 



(4.7) 



Notice that the consistency of the chiral power counting scheme requires that the mesons masses 
appearing in ( [4.6|) be set to their isospin limit; the difference being of higher order, thus, beyond 
the accuracy needed for the present work. 

The isospin breaking term in (|4.6| ) generated by the difference between and is given by 

M^Mk 



-t;+t 



2887r2F2F2 



23 047r2(M|. - M^)^ 



+ 
+ 



4M2. - M2 



Mi-M2 



[36M|. + 563M2M^ - 135M^ - 5767r' M^Mf^B{MK, ±M^)] 



27M^(8Ml + M^)\n 



Ml 



M 2 M2 - 

+ (9Mi - 274M2m|- + 9M^) In ^ + (18M^ + AM^M^ - 9M^) In ^ 



+ 
+ 



M2Mi 



2887r2F2F2. UMj, -M2 

1 r M2 



[68M|: - IQM^ - 127r2(32Mi - 5M^)B{Mk, ±M^)] 



128(M|. + M^) In ^ + (74M| - 43M2) In ^ 1 | . 



M2 



(4.8) 
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The isospin violating term due to the electromagnetic difference between charged and neutral 
mesons masses squared was found to be 



6 



M^Mk 



"11527r2F2F|. 



4715. 



64 A/f2 

^ [1 -187c^ B{Mk,±M^)] 



4M| - M2 



+ 



Mi-M2 



- A71A^kS-^ - 38M^ - 288%^ Mf^B{MK, ±M^) 
Ml . . . , , M2 



+ 27 {Ml + 9M2) In ^ - 2{lQlMl + IIM^) In ^ + (43M| + SIM^) In ^ 1 1 



+ 



MlMl 



78 + 921671^^^ 



11527r2F2F| Ml - M2 
8647r2i3(M;^, ±M^) + 81 In ^ + 110 In 

11^ /i2 



llln^ 

/i2 



M2 



1 



Mi 



11527r2F2F2 Ml- - M2 



64 



M2 



M 



2 ^ 



+ 92167r^[2L^ + - 2(2L^ + L^)] + 63 In ^ - 5 In ^ 



Mi 



64 



5767r2F2F2. l4M2 -M2 



[l + 187r2i3(Mi^,±M^)] 



Mi - M2 



12 



M 



M2 



(4.9) 



+ 92167r2(L"-2L"-Ln + 2887r2S(Mx,±M^)-271n^-341n- , 

with 5ij representing the Kronecker symbol. Finally, virtual photons induce the following 
correction 

, M^K 



6 



2 



1447r2F^F;^ Ml 



l^{6[5Mi + IIM^^ - 3847r2(Mi - M^L;] 



1287r2(Mi - M^){3K{ + - 18K; - 9KI + AK; - 5K;) + 



Ml -Ml 



X 



M2 M2 

(llMi + mlMl + 12M^) In ^ + (4Mi - 39M|.M2 + 3M^) In ^ 



+ 
+ 
+ 



MlMl 

167r2F^Fx Ml 
1 



Ml - M2 



2Mi 



2>F^Fk Ml - Ml 



1^{36 - ?,2^\2Kl - Kl - AK[, - 4K[,) 

M2 M2 " 

(9M|. + 23M2) In ^ - 5(3M| + 4M2) In 
^ [12K; - 6K; + mi + Kl + 77^6' - 12i^8 - 6^[o - ^Ki) 



Ml 



487r2F^Fx Ml - Ml 



12 



- 327r2(12K2^ + Kl + 77^6^ - 127^8" - 18K[o - 18/^[i) - 



9Mi Ml 



Ml - Ml 



.(4.10) 



It is important to stress that the virtual photons contribution, as can be seen from (f4.10|) , is safe 
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from infrared divergencies. The dependence on appears at higher orders of the expansions 



and ([4 .21 ) in power of q. This fact was aheady noticed in ||20| 



For the numerical estimation of the isospin breaking effects we use M^o = 134.976 MeV, M'^o = 
497.672 MeV and the values of the L['s corresponding to set II of Tab. are our preferred ones. 
Concerning the LEG Lg entering the expressions for the electromagnetic form factors of the 
pion and the kaon, its value will be taken as Lg = (6.9 ± 0.7) ■ 10~^. With respect to the EEC's 



in the electromagnetic sector, their central values correspond to the ones quoted in |^ where 
large Nc arguments ||4^ and resonances saturation were applied to pin them down. Moreover, 
an uncertainty of ±1/(167?^), coming from naive dimensional analysis, will be attributed to 
each of the K^^s. In order to compare the order of magnitudes for the scattering lengths and 
the isospin breaking corrections to them at NLO we recall the following combinations from 
Tab. ^ corresponding to the renormalization choice = F^^Fk, 

i (2aJ^ + aj^^ = 0.096 ± 0.012, aj^ = -0.0651 ± 0.0125 . (4.11) 

Having all this, we obtain for the process 7t~K~^ ti~K^ the following estimation for the 
isospin breaking effects 

1.2 10-^ 



327r F^Fk 
1 e 



327r V3 
1 



327r 
1 



327r 

From this we deduce that 



5+-'+- = (2.3±0.1) 10-\ 
2Zoe2F2 5+-/- = (1.4±3.9) 10-\ 
e2 5+-'^" = (-3.8±31) 10-^ 



ao(+-;+-) = 0.097 ± 0.013, (4.12) 

1/2 3/2 

indicating that the isospin breaking effects induce a correction on the combination 2aQ + Oq 
amounting to 1%. As for the process ti^K^ ti^K^ , the results are 



A. 



327r F^Fk 
1 e 



1.2 10-3, 



327r 73 
1 



= (-1.4 ±0.1) 10 



-4 

5 



2Zoe2F2 5++'++ = (-3.3±3.9) 10-^ 



327r --^^^^o-^oe 
1 



327r 



6^(5++'++ = (12.7±31) 10-^ 



and lead to 



ao(++;++) = -0.0631 ± 0.0129. (4.13) 
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Although the size of the correction on is shghtly bigger (~ 3%), it is the former combination, 
(|4.12| ), that is expected to enter the expression for the 2S — 2P energy levels shift for Ktc atoms. 
Before concluding let us comment about the effect induced on these results by deviations from 
the Gell— Mann-Okubo relation and Dashen theorem. If one uses the physical eta mass instead 
of Eq. ( [4. 51 ), the central value in ( ^.12| ) increases by ~ 0.0002. On the other hand, variations 
of 7 in ( ^4.4| ) between 1 to 1.84 cause deviations of order ~ 10^^. Accordingly, deviations from 
both Gell— Mann-Okubo relation and Dashen theorem induce corrections which are beyond the 
accuracy we ask for. 

5 Conclusions 

This work was devoted to the study at NLO of isospin breaking effects on the combination 

1/2 3/2 

2ao + of S-wme ttK scattering lengths which is relevant for the 2S — 2P energy levels 
shift of Kit atoms. We first gave analytic expressions for the scattering lengths which allowed 
to evaluate them using several sets for the values of the LEC's and with both replacements 
Fq — >• Ft^Fk and Fq —> F^. The particularity of the sets in question is that set I and set 
II for instance were determined by fitting the same experimental data to one- and two-loop 
ChPT predictions. We obtained for the value of the afore-mentionned combination a variation 
amounting to ~ 17% due to the choice for the replacement of F^. As for the choice between 
the two sets, it induces ~ 10% variation. These numbers can be used as indicators for the size 
of the NNLO corrections as well as for the rate of convergence of the chiral expansion in powers 
of rris. We pursued by performing, in parallel with what has been done for the tttt scattering 
case, the threshold expansion of the scattering amplitude for the process n^K^ 7i^K~ in 
the presence of photons. After subtracting singular pieces, we evaluated the isospin breaking 
corrections to the combination in question and noticed two features of them: i) in spite of the 
fact that the NLO corrections of both strong and electromagnetic origin are bigger than the 
LO one, they cancel each others at the end. ii) Their central value represents ~ 1% of the NLO 
value for the combination 20^^ + a^^'^ and the uncertainty embarrassing them is negligible with 
respect to the one coming from the LEC's in the strong sector. 
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A Scattering amplitude 

The amplitude ( p2.5| ) will be devided into seven parts depending on the nature of the corre- 
sponding Feynman diagrams 

^+-;+- _ _y\^born _|_ y^tadpole 
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I one- photon _|_ ^^^form factor 



(A.l) 



We distinguish three contributions to the born part of the amphtude (diagram (a) in figure |l|) 



The tree contribution, denoted by A^I?)."™, accounts for the four-meson vertex derived from the 
leading order (LO) lagrangian. Bare masses and mesons wave functions being renormahzed, 
the result is the following 



M 



born 



+ 



+ 
+ 



P2 
-^0 



+ 2e 
4 



6 

2{Ml + 2Ml)Ll + {Ml + Ml)Ll 

-1 Lfin!!!i + ln^ 

87r2 167r2 V My Ml 



2~ 



M|.± r 2 2e . , 16 / e \ 



F2 
-^0 



2e 

{Ml + 2Mi)(2L^ - LI) + Ml{2Ll - L 



V3 

{Ml-Ml){2Ll-Ll) 
1 



{Ml - Ml)Ll 



— e 



-{<oKl + QKl + hKl + hKl - 6K7 - 150K'^ - 2K; - 20Klo - 18K 
9 



Ml. 



1 

3' 



(M^2 + 2M|)(2L^ - LI) + M^(2L^ - U,) 



^-A2F^f^^-^{QKl + QK; 
+ 54K; - 27KI + 5K; + hKl - 6Kj - 78K; - SK; - 134K[o 



126K'^ 



r 3M2 10 



6Ff 
16 



3(M^ + 2Ml)Ll - 3MIl; + 2{Ml + 2M^)Lg + (M^ + MI)L; 



(A.2) 



As usually, the tadpole integrals read 



Hp = Ml Hp 



Ml , Ml 
In — ^ 



327r2F2 



and the difference between the charged and neutral mesons masses is symbolized by 

Ap = Ml± - Mlo . 
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Note that an infrared divergent piece appears in ( [A.2| ). It comes from the charged mesons wave 
functions renormahzation and is regularized by assigning a fictitious mass to the photon. 
The counter-terms contribution A^^'J™ issues from the NLO lagrangian of the strong sector 



and is put in the following form 



Pi 
V4 



1=1 



8(2M2± -t)(2M^± -t) 



4 



U 



2(2M2± - t)(2M|.± - t) + 2{J:^^k+ - s) 



Mi + 14M; 



24e 



(Mi - Ml) 



(2Ml^ -t)- -{12,Ml + 2Ml){2Ml^ - t) 



-{Ml + 2Ml){i:^-K^ 



-{Ml + 2Ml){i:^-K^-u) 



-^(SM^o + Ml^ - A^)(2M|:± -t)- ^(M^o + 3Mi± - 3A^)(2M2± - t) 













8 




3 


V7 


= 0, 




8 




3 


where 





4 

3' 



u 



2e 



2Mi + IhMlMl + M^- -^{2M^ + IIM^M^ - 13M^) 

v3 



(A.3) 



Finally, A^g2p2 represents the counter-terms contribution of 0{e'^p'^). It springs from the NLO 
lagrangian in the electromagnetic sector and reads 



+ 



2e' 
27Fi 

2e^ 
27Fi 

4e^ 
27Fl 

27Fi 



{12KI + 12KI + 54^3^ + 27KI + imi + imi){i:^K - u) 
{mi + QKl + UKl + 27 Kl - AKl + SOKgOl^.i^ - s) 
(3ir[ + b7Kl + 77^5^ + 347^6^(^.7^ - t) 
3{Ml + MDK^ + 3(31M2 + A'iMl)Kl 



+ {ml + MD^g" + (94M2 + 91Mi)i^[o + 90(M2 + Ml)Kl^ 



(A.4) 



The tadpole-type part of the amplitude is furnished by diagram (b) in figure |T|. Its result is 
shown by keeping apart individual contributions from each meson loop 



^yy^ tadpole 



36F2 



6m - t + 6A^-i^+ - 24A^ + ^(3t + 2Ml + 4M|) 

V3 
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+ 



6u-3t- 6M^± + 2M|.± - 20 - ^{t - 2M^ + 4M^) 

v3 



+ -^(12M + 3t-8M^± -12M^± -72A,) 



+ 



2t + 2M|.± + M^i - 3A^ + ^(M| - M^) 

v3 



+ T^(12M + 3t-12M^^±-8M^±-72A,). 



(A.5) 



Concerning the unitary corrections following from diagram (b) and its two crossed in Fig. |I], 
they will be separated according to the channel specifying each crossed diagram. Moreover, 
the contribution of each channel will be divided into parts labelled by the kind of particles 
propagating inside the loop. For instance, in the s-channel 



where 
M 



s-channel 

TT~K + 



\ J s-channel 



\ A s-channel 



\ A s-channel i \ a s-channel i \ a 



^1 6FoV±(s - A^-K+ + 4A^) 



36F, 

+ {s + A^-K++QKfB{Ml±,Ml±-s) 

+ 2{s - ?,K-K+){s + K-K+ + 6A,)Si(M2±, M|:±; s) 

+ {s--iK-K+?B2i{Ml^,Ml^-s) 

+ 2s 



+ 

+ 2 



2Ml± - t + A{u-t-/\^-K+)\ B22{Ml±,Ml±- s) 
^-,[ 2FoVxo [l + ^ (3s - ?,Ml, - Ml 



s - S^o^o + ^(s - ml + IMl) B{Mlo, Ml,- s) 
V3 -I 



s-T.^OKO + -^{s-ml + lMl) X 

s - K~K+ + 4^(^ + 3A^A') 1 5i(M^o, Ml,- s) 
V3 -I 

r e 1 2 - 

+ [s- A,-i^+ + — (s + 3A,A')J 52i(M2o,M|.o;s) 

2 

+ 2s 2Ml±-t + ^{t-2u + 2Ml) B22{Ml, , Ml,-, s) 
L V3 -I 



1 



24Fo^ 



3e \ 



3e 



2FoVko ( 1 - 3s - ?>M^ - Ml, -—{3s -Ml- 3M^ 



V3' 



+ \s + 7M^- ml, - ^ (3s + 1 7 Ml - 23M|) B (M^ , Ml, ; s) 
L v3 -I 



+ 2 



3c 

s + 3A^-A+ - -^(s - A^a) 



X 



+ 



s + - 9M|^o - ^(3s + llMl - 23Ml) BiiM^ Ml,-, s) 
V3 -I 

s + 3A,-A+ - - ^-k) \ B2i{Ml Ml,-, s) 
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+ 2s 



6e 



Au-bt + AMl± - 2M|± - — (t -2u + 2M^) B22{Ml M^o; s) 



(A.6) 



The loop functions figuring in ( |A.6D and in what follows are discussed in Appendix |B], the 
quantity Apg stands for 



ApQ = M| - . 



Similar notations hold for t- and w-channels 



M 



t-channcl 



A ^t-channel i « ^ 



t-channel 



t-channel 



+ M 



M 



M-channcl » /(ii-channel 



t-channel i a ^ t-channel i « ^ t-channel 



Ml 



-K+ 1 



(A.7) 



which individual parts are given by 

1 



A /I tchannel 



36F(f 



6e 

7!^ 



+ 



t-channel 
VV 



72F(f 



12^2 ( 1 



6e 
2e 

7! 



Mi 



+ 



2e 



9t - QM' - 2M^o - ^(9t + 8M^ - 20M^) 5(M^ M'; t) 
V3 -I 



t-channel 



|2FoV(5t-6M2-2M, 



t-channcl 



A ^ t-channel 



+ 2FoV^(5t - 4S^^) + (3t - 4M2)(3t - 3M^ - M^)B{M^, M^; t) 

M^i-l]{s-u) + F^fi^± [5t + 24A^ + 3(s - m) 



167r2 



6 



+ 



^ t(t + - 3{s - u) {mI^ - ] B{Ml^,Ml^-t) , 



+ 



36Fr 
9 
4 



2 



t-channel 
K+K- 



3(s - u) (m^„ - ^) ] 5(M|o, M^o; t) 
^1 FoV± [5t + 3(s - n) + 24A^] 



2 I <± 



(s-m) 



+ 



9 
Li 



t{t + 8A^) - 3(s - ^i) [My - - j J B{Ml^,Ml^- 1) , 



A ^ M-channel 
■'^■K-K+ 



^{ - 6F^fiK±{^n-K+ -u + 2A^) 
9i^0 L 

+ (2M2± +M|.± -m + 3A,)25(M2±,M|^±;m) 
- 2m(2M2± + M|.± - M + 3A,)Si(M2±, M|±; m) 



(A.8) 
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With regard to the one-photon contributions, they will be classified with respect to the topology 
of the exchanged photon 

one-photon .y one-photon , ly one-photon , ly one-photon , »y one-photon 

' 'vertex-leg 's-channel 't-channel ' 'u-channel • 

We distinguish a tadpole-type contribution schematized by diagram (e) in Fig. |l] and for which 
we found 



9p2 

1 /one-photon 

vertex-leg SF^ 



1 



47r2 



1 



47r2 



We also distinguish a unitary contribution given by diagram (d) in Fig. |l|. The expressions 
relative to the different channels defined by the exchanged photon topology read 



« J one-photon 
s-channel 



3F2 I 2 



F^'jl^ Sit - u) - 12{J:^k - u) + s - 2A^K 



+ -F^fiK [3{t-u)- UiE^K -u) + s + 2Ml - QMl 
+ {s- ml - Ml)B{Ml Ml, s) + {?>A^K - s)B^{Ml M^, s) 



+ 2(T.^K - s) 



3{u -t)+ 2i:^K - s G-^is) + 8A^k{^^k - 



M 



one-photon 
t-channel 



^2 



2F^jl^ [6(S,,, -u)- -{2Ml - t) 



- ?,{T.^K-u)B{MlMlt) 



1 



+ 
+ 



167r2 
3 



-9(S^;,-^z)--(t-8M2) 



Q(2Ml - t){^,K - u)GUt) - -is - u)iAMl - 3t)G;,(t)} 



6F,' 



2F^fiK [Q{^.k -u)- -(2M| - t) 



167r2 
3 



1 



-9{^^K - u) - -{t - 8Ml) 



« ^one-photon 
li-channel 



+ 6(2M| - t){J:^K - u)GKK{t) - 2 - - 3t)G^x(^)} 

= -^{P^jl. (5w - 5Mi - 7Ml) + F^jlK i5u - 5Ml - 3Ml) 

+ 2{Ml - u)B{Ml Ml, u) + 2uBi{Ml, Ml, u) 

- 4(S^x - n)(2S^x - u)G:^j^{u) - 4:A^k{^^k - u)G^j^{u) 

+ 6(S^A' - ufG^K^u) + 6(S^A- -u)(^ ) } • 

Finally, the result for diagram (f) in Fig. |l] can be put in the following form 



(A.9) 



M 



form factor ^2 



U — S 
t 



i + r,(t) + rx(t) 



where and F^- represent the electromagnetic vertex functions for the pion and kaon respec- 
tively m 



t 



2L 



, Ml I, Ml 
q t: I In + - In 
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3Fi 



3Fi 



T?2 
^0 



{Ml 



2LI 



~ 4 
1 

96^ 

t. 



B{Ml,Ml, 



1 t 1 
16^ 12. 
1 t 
16^ r2J 



In 



K 



+ - In- 



Ml 

/i2 



{M'^--)B{M'^,M'^,t) + 



{Ml 



)B{Ml,Ml,t) 



1 t 1 
167r2 12. 

1 t 
^ 16^ r2J 



(A.IO) 



B Loop functions 



The loop functions we used in our calculations follow the notations of |^ and which analytic 

2^ . For completeness we recall the expression for the following 



expressions were given m 
function 



B{x,y) 



A/(ar^y)(2aH^ 



{ 



12TT^{x + y) 
arctan 



y/{x - y){2x + y) 



2{x - y) 



+ arctan 



X + 2y 



2^{x-y){2x + y) 



(B.l) 



For the threshold expansion of the processes analyzed in this paper, we needed the following 
expressions for the two-point functions defined in 



ReB{M^i,Mli,Ml 



Re B{Mlo,Mlo,Ml 



1 



167r2 
1 

167r2 



M^±Mk± , 
1 =F TTT^ In 



Ml. 



M\ 



M^±Mk± , 

1 =F TTTT- In 



Ml. 



M\ 



Ml_ 

Ml. 

M^ 
Ml. 



+ 
+ 



A. f 2 

327r2lM|-M2 



2t 



K 



M^ 



{Ml - Miy 

f 2 

327r2 I M| - 



[2(2Mi + Ml) - {Mk ± M^f] In 



2t 



+ 



M2 V" ' Mk 

r2 I o A #2 



[2(M^ + 2M^) - {Mk ± M.)^] In 



Re B{MlMlo,Ml 



{Ml - MIY 
^ [1 + IQtt'^B{Mk±,±M^±)] 



Ml 
Ml 



Ml 
Ml 



+ 



167r2 
1 



487r2 Ml± - 
1 A, 



M2, 



967r2 Ml - M2 



(5M^± T 2M^±){2Mk± ± M^±) In 
2Mk±M^J ^ ^' 



M^ 

Ml. 
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+ 2 



1 



MkTM^ 
2Mk T 



967r2 M| - M2 
V 2MkTM^ 



+ 



Ml 
Ml - Ml 



In^ 

Ml 



(B.2) 



and evaluated at 



M± = M^± ± Mk± . 

Finally, we quote the expression of the infrared divergent three-point function whith the cor- 
responding cut structure needed for the calculation of diagram (d) in Fig. with the various 
channels defined in Eq. ( [A.9| ) 



2Li2 



In 



ApQ + XlQip') 



p2 + ApQ - Xloip^) 



- 2Li, 



p^ 



A 



PQ 



f - ApQ + Xi^ip^) 



In 



[A 



PQ 



'PQy 



[ApQ + XIq{p^W - p^ 



mQ{p 



X 



In 



[p' - ^Up')f - 
[p' + 4q{p')?-aIq^ 



+ 2z7re [p^ - {Mp + Mgf] 



(B.3) 



In this expression, B stands for the Heaviside function, Xpg the Lehmann function 

XpQip') ^ [p' - (Mp + Mq)'] [p' - (Mp - MqY] , 
and the dilogarithm function is defined as 

ln(l - t) 



Li2 (z) = — / dt 



t 



Notice that in the equal mass limit {Mp = Mq), the function (|B.3|) reduces to the one given 
in Pol 
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